HUMAN NEUROSCIENCE 

Mom feels what her child feels: thermal signatures of 
vicarious autonomic response while watching children in a 
stressful situation 

Barbara Manini^-^*^, Daniela Cardone^-^*^ , SJoerdJ. H. Ebisch^-^, Daniela Bafunno^, Tiziana Aureli^ and 
Arcangelo Meria '-^ 

' Infrared Imaging Lab, Institute of Advanced Biomedical Technologies, G. d'Annunzio University, Chieti-Pescara, Italy 
^ Department of Neuroscience and Imaging, G. dAnnunzio University, Chieti-Pescara, Italy 

Maternal attunement with an infant's emotional states is thought to represent a distinctive 
feature of the human primary bond. It implies the mother's ability of empathizing with her 
child in order to fulfil the child's needs in an immediate and appropnate manner. Thus, 
it is particularly involved in stressful situations. By assuming that maternal attunement 
embodies a direct sharing of physiological responses with the child, we compared the 
autonomic response of mothers observing their own distressed child with those of other 
women observing an unknown child involved in an ecological distressful condition (mishap 
paradigm). The hypothesis was that the adult's response was more attuned with the 
child's response in the former group than in the latter group. The autonomic response was 
non-invasively evaluated through the recording of the thermal facial imprints by means of 
thermal infrared (IR) imaging. Nine mother-child dyads and 9 woman-unknown child dyads 
were studied. We found marked similarities between the facial temperature dynamics 
of women and children along the experimental procedure, thus providing evidence for a 
direct emotional sharing within the adult-child dyad. The evidence for common dynamics 
in the time course of the temperatures was assessed through correlation analysis and, 
nevertheless, resulted stronger in the mother-child dyads than in the other women-child 
dyads. In addition, temporal analysis showed a faster response in mothers than in other 
women, thus confirming our study hypothesis. Besides confirming the extraordinary 
capability of IR imaging to preserve ecological context in the study of social or non-verbal 
interactions, these results suggest that maternity appears to potentiate the emotional 
attunement with the child. Although based on preliminary results, this study opens new 
perspectives in the study of the factors modulating vicarious socio-emotional processes. 

Keywords: emotion, vicarious responses, emotional shiaring, motlier-child synclirony, motlier-chiild bond, 
autonomous nervous system, IR thermal imaging 
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INTRODUCTION 

Maternal attunement with a child is a key element of secure 
attachment relationships (Bowlby, 1958). A mother who is capa- 
ble to share affect with her child, to empathize with his/her 
emotional needs and to appropriately respond to his/her requests, 
allows the infant to perceive a sense of being accepted and recog- 
nized, which facilitates social adjustment and a positive psycho- 
logical functioning (Bowlby, 1958). Conversely, the lack of such 
a sensitivity could lead to less favorable outcomes in the child's 
subsequent development (Sullivan et al., 2011). Although there is 
a general agreement on the importance of maternal attunement, 
little is still known about the biological basis of that ability. 

Maternal attunement with offspring is especially involved 
when mothers are facing infants' distress. As shown by a large 
amount of empirical evidence, mother's alertness and arousal 
increase because of the baby's distress signals (Swain et al, 2007), 
allowing to immediately provide the help needed for the infant's 
recovery (Bell and Ainsworth, 1972; Frodi and Lamb, 1978; 



Mills-Koonce et al., 2007). As found by neuroimaging research, 
watching a child's emotional expressions specifically activates 
some brain regions, like the anterior insula, the amygdala and 
the mirror neuron system (Rizzolatti et al., 2001; Gallese et al., 
2004). This activation occurs mostly when a mother observes her 
own child rather than an unknown child (lacoboni et al., 1999; 
Lenzi et al., 2008). Since the anterior insula is considered the relay 
between the action representation (mirror neuron system) and 
the emotional processing (limbic system) (Lenzi et al, 2008), the 
above activation could relate to the mother understands of her 
own child's distress signal, probably to the aim of responding suc- 
cessfully and appropriately. Particularly, the infant's crying — the 
most powerful distress signal — activates the same brain regions 
involved in attention, emotional attachment and in the process of 
integrating autonomic states with social behavior (Seifritz et al., 
2003; Swain et al, 2007, 2008; Kim et al, 2010; Swain, 2011; 
Swain et al., 2011). The infant's distress also elicits autonomic 
responses in the own mother. In particular, blood pressure as 
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index of aversive stimuli, heart rate as attentional and defensive 
response, and skin conductance as index of arousal, increase in 
mothers witnessing their infant's crying (Frodi and Lamb, 1978). 

Beyond the ability to provide the right intervention in distress- 
ing situations, synchronization of the mother's responses to the 
infant's signals in typical dyadic interactions can be considered 
a key aspect of sensitive parenting, as it implies the prompt- 
ness of the mother response (Bell and Ainsworth, 1972) and the 
adaptation moment by moment to the child's emotional states 
(Noriuchi et al, 2008; Guedeney et al, 2011). Repeated expe- 
rience of well-synchronized and appropriate interactions allow 
mother and child to become sensitive to each other's physiolog- 
ical and behavioral cues (Fleming et al, 1999; Feldman, 2007; 
Feldman et al, 2011, 2012) and capable to perceive the other's 
behaviors as a response to their own behavior, which contributes 
to the formation of a unique bond between them (Stern, 1985; 
Mogi et al, 2011). 

The autonomic nervous system seems to represent an ele- 
mentary mechanism supporting emotional synchrony between 
mother and infant. According to Forges (1998, 2003a,b), the 
emergence of parenting behavior in mammals is linked to the 
development of the Polyvagal system (Feldman and Eidelman, 
2007). In his opinion and with reference to the human realm, 
the decrease of heart-rate variability due to vagal-tone suppres- 
sion in response to a stressful event (Forges, 2003a,b) would 
facilitate complex behaviors such as attention, orientation and 
the maintenance of calm states, which are required for early 
formation of secure social bonds and for achieving the more 
advanced coordination of social signals that underlie human 
interpersonal interactions. Supporting that hypothesis, research 
on early infancy found that the mother's and her child's vagal 
tone are inter-related (Feldman et al., 2010), thus suggesting that 
the degree of parasympathetic control during social engagement 
is shaped by co-regulatory processes. As suggested by Feldman 
(2007); Feldman et al. (2012), biological synchrony in inter- 
personal interactions shows online sensitivity to the partners' 
ongoing behavior. 

The physiological side of the mother-infant bond has been typ- 
ically observed in stressful situations. Mother's and infant's Corti- 
sol reactivity were found to be strictly associated in the stiU-face 
paradigm, an interactive situation eliciting some degree of stress 
in infants (Haley and Stansbury, 2003; Feldman et al., 2010), 
and while playing a challenging- for-the-child game (Sethre- 
Hofstad et al, 2002). All together, these findings provide pro- 
totypical examples of the bio-behavioral synchrony underlying 
interactions, as suggested by Feldman (2007). Accordingly, an 
autonomic-visceral synchrony between mother and child was 
found in our previous study (Ebisch et al., 2012), showing a 
significant parallelism between mothers' and children's facial 
temperature variations when mothers observed their own chil- 
dren involved in a stressful situation. Since skin temperature 
is mediated by the autonomous nervous system and varies in 
response to emotional stimuli originated from external environ- 
ment (Anbar, 2002; Merla and Romani, 2007; Nakanishi and 
Imai-Matsumura, 2008; Shastri et al, 2009; Nhan and Chau, 
2010; Kuraoka and Nakamura, 2011), that study showed for the 
first time a direct sharing at the autonomic level between mother 
and child associated with an experience of affective attunement. 



The present study aimed to deepen the results found by our 
previous one (Ebisch et al., 2012). In that study we demonstrated 
autonomic thermal synchrony in mother-child dyads when the 
mother watched her own child experiencing a stressful situation 
elicited by the mishap paradigm (Cole et al, 1992). Thermal 
autonomic responses were assessed through the use of high- 
resolution thermal infrared (IR) imaging (Merla and Romani, 
2007). In the present study, using the same paradigm of the 
previous study, we compared the thermal autonomic responses 
between two groups of dyads: mother-her own chUd dyads and 
women-unknown chUd. Since, according to previous research 
(Seifritz et al, 2003; Swain et al, 2007, 2008), the experience of 
one's own child distress seems to influence the women's reac- 
tions to the infant's emotional needs, we expected to find different 
autonomic responses, on either the intensity or the time scale, 
across the groups. Specifically, we hypothesized that the former 
group of dyads would attune more than the second group, thus 
showing a stronger shared response. Alternatively, the two groups 
could exhibit the same physiological response, thus suggesting 
that the such a response to the child's distress is a general reac- 
tion, independently of the specific bond with the distressed child. 
Moreover, according to evidence for the contingency of mater- 
nal responses as one of the specific characteristics of maternal 
sensitivity (Feldman, 2007; Feldman et al., 2012), we expected 
to find more ready and prompt responses in the first group of 
dyads. 

METHODS AND MATERIALS 
PARTICIPANTS 

Fourteen children (seven male, age 39-45 months) participated 
in the study. All of the children were born at term and had a 
typical psychological and physical development. The adult sam- 
ple was composed of eighteen women divided in two groups: 
the first one included nine mothers (age 25-43) who watched 
their own child during the experiment ("mothers," M); the second 
sample was composed of 9 women (age 23-38) ("other woman," 
OW). Four of these had preschool children not participating to 
the study, while the other five were not mother. The data about 
six out of nine mother-child dyads were included in our previ- 
ous study (Ebisch et al., 2012). The two groups of women were 
matched for socio-economical status and study degree. Inclusion 
criterion for adult subjects was the absence of any overt phys- 
ical, psychiatric or psychological disease. All participants were 
asked to refrain from heavy physical activities and consumption 
of coffee, cigarettes and vasoactive substances for at least 2 h prior 
to the measurements, and to avoid cosmetics on their faces at 
the moment of the experiment. The study was approved by the 
Local Ethics Committee. Written informed consent was obtained 
from all the participants in agreement with the Declaration of 
Helsinki. 

PROCEDURE 

Prior to testing, each subject was left to acclimatize for 10-20 min 
to the experimental room and to allow the neutral condition skin 
temperature to stabilize. The recording rooms were set at stan- 
dardized temperature (23°C), humidity (50-60%) without direct 
ventilation. The subjects sat comfortably on a chair during both 
the acclimatization and the measurement periods, without any 
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restriction of body movements. Before the start of the experiment, 
the children underwent an adequate familiarization period to ease 
the psychological habituation to the setting and the experimenter, 
first in presence of their mothers, then followed by neutral inter- 
action with the experimenter alone. During neutral interaction 
between the experimenter and the child, some toys were presented 
to allow the child to feel at ease and to get used to play with the 
experimenter. 

After the neutral interaction with the experimenter, the chil- 
dren were exposed to a potential stressful experience, elicited 
by the "Mishap Paradigm" (Cole et al, 1992). More specifically, 
children were invited to play with a toy, which was previously 
manipulated to break in the child's hands when playing with it, 
thus suggesting that the child accidentally broke the toy. The toy 
was introduced by the experimenter as her own favorite. Distinct 
phases could be distinguished in the paradigm: (1) "presenta- 
tion" (the experimenter demonstrated the toy); (2) "playing" (the 
chUd played with the toy, while the experimenter left the room for 
1 min); (3) "mishap" (child "broke" the toy); (4) "re-entry" of the 
experimenter (the experimenter did not say anything for 30 s and 
merely looked at the broken toy); (5) "soothing" of the child (the 
experimenter cheerfully indicated that the toy could be fixed and 
that the breaking was not the child's fault). 

In order to perform the analyses, the above mentioned phases 
were grouped in conditions. The neutral interaction defined the 
"neutral condition"; the presentation and playing phases together 
formed the "intermediate condition"; the mishap, re-entry and 
soothing phases together formed the "experimental condition." 

The "mothers" and the "other women," naive about the spe- 
cific content of the experiment, were invited to silently observe 
the children-experimenter interaction through a one-way mirror 
from a separated room. It was possible that two women (one for 
each group of dyads) watched at the same child together, though 
they could not see each other or interact as a screen was placed 
between them. 



Facial thermal images for all of the subjects were recorded 
along the whole experimental procedure (acclimatization, neu- 
tral, intermediate, and experimental conditions). 

MATERIALS 

Thermal IR imaging was performed by means of three digital 
thermal cameras FLIR SC660 (640 x 480 bolometer FPA, sensi- 
tivity: <30mK @ 30°C). The acquisition frame rate was set to 
15 Hz for each thermal camera. 

Two remote-controlled video-cameras (Canon Vc-C50iR) 
were used to film the child for the behavioral analysis. Video- 
signals were sent to two video-recorders (BR- JVC) and mixed by 
a Pinnacle system (Liquid 6). Subsequently, the movies were pro- 
cessed through a specialized software (Interact Plus, Mangold) 
that allows to code behavior in synchrony with the ongoing 
movies of the children during the experiment. The toy pre- 
sented to the children in the "Mishap Paradigm" was a black 
and white robot with a height of approximately 20 cm. When 
turning on the robot with a switch on its back, it started to 
walk and play music. Both hands of the robot could be opened 
and closed by means of pressing/relieving a button. One of the 
hands of the robot were prepared to break when manipulated 
by the child. The robot could be repaired only by the experi- 
menter. The toys presented during neutral interaction between 
the experimenter and the child were a puzzle, a magic wand 
and 3-D book. 

DATA ANALYSIS 
BEHAVIORAL DATA ANALYSIS 

Following from the notion that different combinations of rel- 
evant signs in a mishap situation may be indicative of guilt 
(Barrett et al, 1993; Kochanska et al, 1995), we coded the pres- 
ence or absence of the child's reactions into five categories (see 
Table 1): gaze and eye, bodily tension, arms, repair, and verbal- 
izations (Barrett et al, 1993; Kochanska et al, 2002; Mills, 2003; 



Table 1 | Expressive features coded for the scoring of child reactions. 



Category 


Features 


Description 


Gaze and face''''' 


Gaze aversion 


The child stares into space, or toward the oblique low, or toward another insignificant object 
(excluding the broken toy and the experimenter). 




Lip rolled-in 


Lower lip rolled-in; corners of mouth drawn. 


Bodily tension'' 


Bodily avoidance 


The child backs up while looking at the experimenter; or moves away from her, toward 
insignificant object, after focusing on her. 




Hunched shoulders 


Relaxed or hunched shoulders. 




Head lowered 


Head hanging or tilted forward. 


Arms'' 


Arms across body 


Arms across the midline, held close to the body (e.g., hugging the body). 




Covering, touching face 


The child covers or touches all or part of the face. 




Fingers in mouth 


Putting a finger or fingers in mouth. 


Repair" 


Trying to repair the object 


The child tries to repair, to fix the toy or he/she manipulates it. It is not coded as repair if the 
action is meant to play. 


Verbalizations''-'' 


Confession 


The child admits to have broken the toy e.g., saying "1 broke it" or "1 pulled this piece off." 




Negative self-evaluation 


The child judges him/herself negatively, e.g., saying "1 am not able to play" or "1 can't do it." 



''Barrett et al., 1993; Barrett, 2005. 
1^ Mills, 2003. 

'^Tangney and Fischer, 1995; Stipek, 1995; Lewis et al., 1998: Kochanska et al., 2002. 
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Barrett, 2005). The duration of the phases was similar among chil- 
dren [M = 35.86 (SD = 17.54) seconds in playing; M = 62.21 
{SD = 14.52) seconds in mishap; M = 48.86 (SD = 16.44) sec- 
onds in re-entry] , and there were no cases of distress so important 
to reduce the duration of the phase in a relevant way for the 
appearance of the behaviors to code. Thus, in playing, mishap 
and re-entry of experimenter, overall distress response was coded 
on a 4-point scale [(1) child was not affected in any way; (2) 
child appears middle distressed as evidenced by one or two sig- 
nal coded; (3) child appear distressed as evidenced by three or 
four signal coded; (4) child is strongly distressed as evidenced 
by five signal coded]. Behaviors were coded with the software 
INTERACT 8. 

Comparisons were made between the playing, mishap and re- 
entry phases, because in these phases the children were engaged 
with the same object and that made more reliable the coding of 
their behaviors. 

Reliability 

Two observers coded all the children for each phase and Kappa's 
ranged from 0.72 to 0.80 (allp < 0.05). 

THERMAL DATA ANALYSIS 

A visual inspection of the changes in facial thermal imprints in all 
subjects was performed to qualitatively investigate the autonomic 
responses of the women and the children throughout the exper- 
iment. This analysis was followed by a quantitative estimation 
of temperature variations of the nasal tip. Referring to the liter- 
ature about thermal signature of emotions (Shastri et al., 2009; 
Nhan and Chau, 2010; Kuraoka and Nakamura, 2011), the nasal 
tip has been proved to be associated with the activation of the 
sympathetic nervous system by emotional and distressing stimuli 
(Roddie, 1963; Merla and Romani, 2007; Nakanishi and Imai- 
Matsumura, 2008; Shastri et al, 2009; Nhan and Chau, 2010; 
Kuraoka and Nakamura, 2011; Ebisch et al., 2012). More pre- 
cisely, temperature decrease of the nasal tip is related to sympa- 
thetic alpha-adrenergic vasoconstrictor activity and the reduction 
of this vasomotor effect is related to the thermal recover due to the 
vasodilatation effect. 

Thermal signals have been extracted through the use of track- 
ing software, developed with homemade Matlab algorithms (The 
Mathworks Inc., Natick, MA). The tracking algorithm is based on 
the 2-D cross-correlation between a template region, chosen by 
the user on the initial frame, and a similar ROI in a wider search- 
ing region, expected to contain the desired template in each of 
the following frames (Tangherlini et al, 2006). In this way it is 
possible to automatically extract the thermal signals in defined 
regions of interest during the whole experiment. The extracted 
thermal data have been filtered subsequently with a low-pass fil- 
ter (fcut-off = 0.2 Hz), to eliminate breathing effects (Ebisch et al., 
2012). 

A comparison was made between the neutral, the interme- 
diate, and the experimental conditions. Because the mishap 
paradigm is an ecological situation the duration of the conditions 
depended on the children spontaneous behaviors. Therefore, the 
timing of the experiment presented some variations within sub- 
ject. The mean duration of neutral condition was 134.62 {SD = 



80.21) s. The mean duration of the intermediate condition was 
139.104 {SD= 115.097) s. For the experimental condition the 
mean duration was 222.327 {SD = 55.79) s. 

STATISTICAL DATA ANALYSIS: INDIVIDUAL DYADS 

The thermal signals of the nasal tip for all the subjects were trans- 
formed in z-scores. A repeated measures (3x3) ANOVA was 
performed on the mean nasal tip temperature of the three con- 
ditions for all subjects. The conditions (neutral, intermediate, 
and experimental) were set as within-subject factor. The groups 
(children, mother, other women) were set as between-subject 
factor. Pairwise multiple comparisons between conditions were 
adjusted for multiple comparisons by means of Bonferroni cor- 
rection. The goal of this analysis was to control the autonomic 
response of the subjects during the entire paradigm. In particular 
our hypothesis was that there were no differences in autonomic 
response of the subjects between neutral and intermediate con- 
dition. On the contrary, it was hypothesized that the autonomic 
response in the experimental condition was significantly dif- 
ferent, compared with both the neutral and the intermediate 
condition. Besides, we wanted to investigate if this autonomic 
response changed according to the specific group of the subjects. 
Pearson correlation coefficients among the time courses of the 
temperature signals were calculated for each women-child dyad 
and for the neutral and the experimental condition. Correlation 
coefficients were then transformed according to Fisher's distribu- 
tion (r to z transform) and a f-test was performed to evaluate 
whether significant differences could be found between condi- 
tions (neutral and experimental) and groups of dyads in the two 
conditions. 

A cross-correlation like analysis on each dyad's temperature 
signals during the experimental condition was also performed to 
test which delay between the signals corresponded their highest 
correlation value. For this purpose, Pearson correlation coeffi- 
cients were computed by delaying the woman's signal with respect 
to the child's signal with 10, 20, 30, 40 s. For each group, the mean 
values of the delays corresponding to the highest dyad correlation 
were computed. A Mann- Whitney [7-test was performed to assess 
possible differences between the two groups of dyads regarding 
the distributions of time delays maximizing the correlations. 

As the "other women" group included four mothers, explo- 
rative descriptive analyses on the "other women" subgroups 
were performed to preliminarily investigate possible differences 
between other mothers and non-mothers. Pearson correlation 
coefficients were calculated on the individual dyads level for the 
other mothers and non-mothers subgroups. Furthermore, the 
temporal modulation of the signal was explored for the two 
subgroups. 

GROUP ANALYSIS 

Since the experiment duration depended on the child's behavior, 
signal resampling was performed to compare the three conditions 
and the three groups (children, mother, other women). Each tem- 
perature signal was resampled in order to obtain 300 data points 
for each condition, equally spared across the time duration of 
the condition (for the experimental condition, 120 samples were 
included for the mishap phase, 60 for the re-entry phase and 
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120 for the soothing phase). The resampled z-score temperatures 
where then averaged for the groups, thus providing an average 
time course for each group and condition. 

RESULTS 

BEHAVIORAL RESULTS 

Behavioral analyses showed that all children experienced an 
increase of distress in response to the mishap. Specifically, the 
repeated-measures ANOVA indicated that the children reacted 
in a different way in the playing, mishap and re-entry of exper- 
imenter phase [_F(i^ 13) = 10.95; p < 0.05]. There was an evident 
increase of distress signs in both the mishap and experimenter 
re-entry phase, compared with playing (f = 4.37, p < 0.01; and 
t = 3.31, p< 0.01, respectively). No significant evidence of dif- 
ference in distress signs was found comparing the mishap and the 
experimenter re-entry phases. Table 2 reports the results of the 
behavioral analysis. 

VISUAL ANALYSIS OF FACIAL THERMAL IMPRINTS 

A visual inspection of the facial thermal imprints was performed 
to investigate the presence of appreciable signs of autonomic 
responses in the mothers, the other women and the children 
throughout the experiment. The qualitative results resembled 
those already found in our previous work (Ebisch et al, 2012). 

A representative example of the facial thermal imprints of one 
mother-other woman-child triad is shown in Figure 1. As to the 
child, no appreciable modulations were detected regarding the 
facial skin temperature distribution between the neutral condi- 
tion and the intermediate condition. However, after the mishap 
a sympathetic reaction could be observed, reflected by a sud- 
den and wide-spread decrease of face temperature, especially in 
the nasal tip as previously found in human as well as macaques 
(Kistler et al, 1998; Merla and Romani, 2007; Nakanishi and 
Imai-Matsumura, 2008; Shastri et al, 2009; Nhan and Chau, 
2010; Kuraoka and Nakamura, 2011; Ebisch et al, 2012). The 
decreased skin temperature in the nasal tip likely reflects periph- 
eral vasoconstriction due to the alpha-adrenergic activity. These 
sympathetic responses gradually decrease after the re-entry of the 
experimenter. During the soothing phase, the temperature of the 
nasal tip soon increased, likely reflecting a withdrawal of the sym- 
pathetic alpha-adrenergic vasoconstrictor effect. Moreover, the 
nasal tip temperature presented an over-response, compared to 
the neutral condition value. 



Concerning the mother, no appreciable modulation of skin 
temperature distribution was detected during the intermedi- 
ate condition. After the mishap as well as after the re-entry 
of the experimenter, the same thermal variations observed in 
the child could be appreciated in the mother. During the 
soothing phase, the mother showed a gradual and general- 
ized increase of facial temperature, re-establishing the neutral 
condition state. Moreover, like the child, she showed an over 
response of the nasal tip temperature, compared to the start 
of the experiment. 

Finally, the other woman did not show any signal modulation 
between neutral condition, intermediate condition, and mishap. 
During the experimenter re-entry, there was a cooling of the nasal 
tip, however, milder (using the same temperature scale and image 
contrast) than that showed by the child and the mother. In the 
soothing phase the other woman showed a recovery of the nasal 
tip temperature to the initial neutral condition value. 

COMPARISONS OF AUTONOMIC RESPONSES BETWEEN GROUPS AND 
CONDITIONS 

An ANOVA (3 x 3) was performed on the re-sampled z-score 
mean temperatures of the nasal tip. A significant within- 
subject effect was found for condition [_F(2, 29) = 5.70;p < 0.01]. 
Bonferroni adjusted pairwise comparisons on the within-subject 
factor "condition" showed no significant difference between neu- 
tral and intermediate condition {p = 1). There was a significant 
difference between the neutral and the experimental condition 
ip < 0.05) and between the intermediate and the experimen- 
tal condition (p < 0.05). The condition x group effect was not 
significant [F(2, 29) = 0.58; p = 0.68]. These results suggest that 
the experimental condition, but not the intermediate condi- 
tion induced a significant modulation on autonomic response 
in all three groups (i.e., children, mothers and other women), 
compared with the neutral condition. 

CORRELATION ANALYSIS 

Table 3 shows the Pearson correlation coefficients calculated for 
all the dyads in the neutral condition and in the experimental 
condition. As we were interested in studying the attunement of 
the autonomic responses of the mothers with those evoked by the 
distressful condition in the child, the non-distressftil intermediate 
condition was not included in the correlation analysis. Pearson 
correlation coefficients were found to be statistically significant 



Table 2 | Means and standard deviation of the categories in all phases coded. 



Categories 



N 



Mean 



Playing 



SD 



Phases 



Mishap 



Mean 



SD 



Re-entry 



Mean 



SD 



Gaze and face 
Bodily tension 
Arms 
Repair 

Verbalizations 
Total scores 



14 
14 
14 
14 
14 
14 



0.57 

0.14 

0.21 

0 

0 

1.71 



0.51 

0.36 

0.42 

0 

0 

0.47 



0.86 
0.36 
0.29 
0.64 
0.21 
2.42 



0.36 
0.50 
0.47 
0.50 
0.43 
0.51 



0.71 
0.50 
0.29 
0.43 
0.21 
2.28 



0.47 
0.52 
0.47 
0.51 
0.43 
0.61 
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FIGURE 1 I Facial thermal imprints of one of the mother-other woman-child triads. 



Table 3 | Pearson correlation coefficient for all the dyads. 



Dyad 


Type of 


Neutral 


Experimental 




dyad 


condition 


condition 


1 


M-C 


0.53* 


0.67* 


2 


M-C 


0.51* 


0.76* 


3 


M-C 


-0.51* 


0.71* 


4 


M-C 


0.49* 


0.84* 


5 


M-C 


0.86* 


0.91* 


6 


M-C 


-0.64* 


0.94* 


7 


M-C 


-0.61* 


0.28* 


8 


M-C 


0.86* 


0.97* 


9 


M-C 


0.06* 


0.23* 


10 


OW-C 


-0.12* 


0.24* 


11 


OW-C 


0.26* 


0.58* 


12 


OW-C 


-0.77* 


-0.03* 


13 


OW-C 


-0.85* 


0.82* 


14 


OW-C 


0.35* 


-0.38* 


15 


OW-C 


-0.24* 


0.75* 


16 


OW-C 


0.29* 


0.40* 


17 


OW-C 


0.93* 


0.71* 


18 


OW-C 


-0.32* 


0.54* 



In the column type of dyad M-C, mother-child dyads; OW-C, other women-child 
dyads. *p < 0.001. 



(allp < 0.001). With respect to the mother-child dyads group, the 
mean coefficient value resulted rmean = 0.17 (SD = 0.60) during 
the neutral condition and during the experimental condition it 
increased up to a mean value of rmean = 0-70 (SD = 0.27). The 



t-test on the Fisher-transformed r-values showed a significant 
difference in the correlation coefficient between neutral condition 
and experimental condition in the mothers-child dyads group 
(f = -3.32, d/ = 8,p < 0.05). 

With respect to the other woman chUd-dyads, in 5 out of 
9 cases there was an increase of correlation in the experi- 
mental condition, compared with the neutral condition. The 
mean value of the neutral condition correlation was Tmean = 
—0.05 (SD = 0.57), while the mean value in the experi- 
mental condition was r^em = 0.40 (SD = 0.39). The f-test 
showed no differences between neutral condition and exper- 
imental condition in this second group (f=— 1.64, df = 8, 
p = 0.14). 

The comparison between the two groups suggested that there 
was no significant difference regarding correlation coefficients in 
the neutral condition (f = 0.76, df = 16, p = 0.45). Instead, the 
t-test yielded a significant difference between the two groups in 
the experimental condition (t = 2.19, df = 16, p < 0.05), reflect- 
ing a higher mean correlation in the mother group, compared 
with the other women group. 

Finally, an explorative descriptive analysis was performed in 
order to assess the differences in correlation coefficient between 
other mothers-child dyads and non-mother-child dyads belong- 
ing to the other women group. The mean correlation in the 
other mother-child dyads was rmean = 0.16 (SD = 0.58) for 
the neutral condition and rmean = 0.60 (SD = 0.58) for the 
experimental condition. In the non-mother-child dyads the 
mean correlation is r^ean = 0.16 (SD = 0.58) for the neutral 
condition and rmean = 0.25 (SD = 0.47) for the experimental 
condition. 
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TEMPORAL ANALYSIS OF THE RESPONSES 

Table 4 shows the Pearson correlation coefficients for the thermal 
signal of the nose tip of adult and child for the experimental con- 
dition, in both the groups of dyads. In order to evaluate whether 
the maximum value of the correlation varied along the time 
course of the experimental phase, i.e., to evaluate possible delayed 
responses, and whether such eventual delays differed between 
the two groups of dyads, we computed the Pearson correlation 
coefficients by delaying the woman's signal with respect to the 
child's one with 10, 20, 30, 40 s, like in a cross-correlation analysis. 
The largest delay with respect to which we report the cross- 
correlation results is 40 s, since a higher delay did not maximize 
the correlation value in any dyads (see Table 4). 

The average value of the maximum correlation coefficients 
for the mothers-child dyad group was rmean = 0.81 (SD = 0.21), 
corresponding to an average delay of 7.8 (SD = 13.94) s. 

One mother-child dyad (dyad #6) reached its maximum corre- 
lation value at a 40-s delay, even though its correlation coefficient 
was already larger than 0.67 and significant at 0 time delay. By 
excluding this dyad from the group, the average value of the 
highest correlation coefficients did not change, but there was a 
significant decrease of the mean delay in which r-value is max- 
imized in the mother-child dyads. Excluding mother-child dyad 
number six reduced the average delay to 3.75 s (SD = 7.44). 

The average value of the maximum correlation coefficients 
for the other woman-child dyad group was rmean = 0.57 (SD = 
0.42), corresponding to an average delay of 26.7 (SD = 20.0) s. 
The Mann- Whitney l/-test on the delay levels showed a signifi- 
cant difference between the two groups (11= 14, p < 0.05), even 
when taking into account the dyad #6 in the mother-child dyad 
group. 



Although aware of the small size of the sample, an explorative 
analysis was performed on the other women group, by dividing it 
in two subgroups: other mothers (i.e., mothers watching not their 
own child) and non-mothers. The mean of the maximum corre- 
lation coefficients in other mother-child dyads was rmean = 0.64 
(SD = 0.21), whereas the mean of the delays was 20 s (SD = 
23.09). In the non-mother-child dyads the mean of maximum 
correlation wa,s t'lnean — 0.52 {SD — 0.56) with ail avcrcigc delay 
of32 s (SD = 17.88). 

GROUP ANALYSIS OF CORRELATION 

Figure 2 shows the groups' average signal. The correlation coef- 
ficients were calculated among the three average signals, for 
the their whole time course and for each condition. For the 
whole time course signal, the correlation coefficient value for the 
mother-child dyads was r = 0.88 (p < 0.001), while it resulted 
r = 0.84 (p < 0.001) for the other women-child dyads. The 
Pearson coefficient values did not differ between the two groups 
in the experimental condition (mother — child: r = 0.89 p < 
0.001; other women— chUd: r = 0.90 p < 0.001). 

DISCUSSION 

The present study focused on the physiological expression of 
the emotional attunement of mothers to their own children. 
We compared the autonomic response, measured through the 
nose tip temperature, of a group of mothers while watching 
their own child engaged in a stressful situation with that of 
another group of women observing unknown children in the 
same situation. According to the exiting literature, reporting 
increased behavioral responses and higher neural activation in 
mothers watching their own distressed chUd, we expected that 



Table 4 | The table shows the Pearson correlation coefficient for all dyads for synchronous and shifted signals. 



Delay 



No. 


Type of dyads 


No delay 


IDs 


20 s 


30 s 


40 s 


1 


M-C 


0.24* 


0.96* 


0.94* 


0.92* 


0.90* 


2 


M-C 


0.28* 


0.24* 


0.21* 


0.19* 


0.16* 


3 


M-C 


0.97* 


0.34* 


0.41* 


0.5* 


0.58* 


4 


M-C 


0.82* 


-0.34* 


-0.35* 


-0.39* 


-0.29* 


5 


M-C 


0.76* 


0.72* 


0.94* 


0.47* 


0.92* 


6 


M-C 


0.67* 


0.76* 


0.78* 


0.83* 


0.85* 


7 


M-C 


0.71* 


0.63* 


0.44* 


0.19* 


-0.03* 


8 


M-C 


0.84* 


0.74* 


0.50* 


0.13* 


-0.20* 


9 


M-C 


0.91* 


0.89* 


0.69* 


0.64* 


0.69* 


10 


OW-C 


0.75* 


0.79* 


0.84* 


0.88* 


0.91* 


11 


OW-C 


0.40* 


0.30* 


0.28* 


0.14* 


-0.17* 


12 


OW-C 


0.71* 


0.67 


0.68* 


0.62* 


0.35* 


13 


OW-C 


0.51* 


0.28* 


0.35* 


0.45* 


0.55* 


14 


OW-C 


0.24* 


0.60* 


0.68* 


0.71* 


0.74* 


15 


OW-C 


0.58* 


0.51* 


0.49* 


0.49* 


0.47* 


16 


OW-C 


0.28* 


0.48* 


0.60* 


0.68* 


0.74* 


17 


OW-C 


0.55* 


-0.52* 


-0.53* 


-0.52* 


-0.46* 


18 


OW-C 


0.82* 


0.79* 


0.86* 


0.88* 


0.98* 



In the column type of dyad M-C, mother-child dyads; OW-C, other women-child dyads. *p < 0.007. Bold value indicates maximum correlation coefficient for the dyad. 
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#Frames# 



FIGURE 2 I Mean signals of mothers, other women and children for all the conditions. 



the two groups would have presented differences with respect 
to either the intensity and/or the timing of their physiological 
reactions. Our results indicate that, for both groups of women, 
the considered autonomic response significantly changed dur- 
ing the experimental condition with respect to the neutral one. 
The emotional attunement resulted slightly stronger and much 
faster in the mother-child group than in the other women-child 
group. 

Specifically, according to the hypothesis that the emotional 
attunement with another person embodies a direct sharing of 
visceral-autonomic responses (Konvalinka et al, 2011), we found 
that the time course of the nose tip temperature confirmed that 
adult women shared such an autonomic response with child 
exposed to a distressful situation, independently of parenting 
relations with the child. When analyzing the single adult-child 
dyad response, the two groups of dyads differed. The autonomic 
responses stronger correlated during the experimental condition 
with respect to the neutral condition for all of the mothers-child 
dyads, while the same did not always happen in the other women 
group. 

Moreover, the highest cross-correlation coefficient values 
between the time courses of the adult-child dyad signals during 
the experimental condition were found at a very short or null 
delay for the mother-child group, whereas a much longer delay 
characterized the other women-child group. 

Group average correlation coefficient values did not signif- 
icantly differ between the two groups. This result may rise an 
apparent discrepancy between the results found at a group level 
and at a single dyad level results. Such a mismatch likely depends 
on the procedure used to perform the group analysis, which 



emphasized the shared signals modulations within the dyads. As 
the experiment was designed to maximally preserve the natural 
and ecological context, there were no a priori limits or timing for 
the time duration of the phases' time duration. Therefore, it was 
not possible to adopt a unique common time frame the whole at 
the group level. 

The last finding suggests that the emotional promptness to 
one's own child's distress is a distinct feature of maternal attitude, 
thus confirming at a physiological level the largely confirmed 
results found at a behavioral level (Stern, 1974, 1999; Brazelton 
et al, 1974; Beebe and Lachman, 1988; Tronick, 1989; Trevarthen, 
1993; Fleming et al, 1999; Trevarthen and Aitken, 2001; Feldman, 
2007; Feldman and Eidelman, 2007; Papousek, 2007; Feldman 
et al, 2012). Since such a faster promptness for emotional vicar- 
ious processes seems to characterize maternal responses, our 
results further suggest to consider it as a basic mechanism for 
organizing early mother-child interactions. The faster autonomic 
response of the mothers supports the idea that the Autonomic 
Nervous System plays a fundamental warning role for the mother 
with respect to the child's emotional shift (Frodi and Lamb, 1978; 
Gunnar and Donzella, 2002). This signal would modulate an 
embodied and pre-reflective sensitivity that helps the mother to 
immediately recognize any shift in the child's emotional needs, 
as well as to promptly soothe the child when distressed (Mills- 
Koonce et al, 2007). 

The explorative, qualitative analysis within the other women 
group suggested interesting preliminary information about pos- 
sible differences between mothers and not mothers looking 
unknown children. The small size of the sample does not allow 
drawing any conclusion, but our results suggest to investigate 
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further whether the other mothers may have a stronger 
autonomic attunement with child and a faster response than 
non-mothers. Such a result would be in accordance with pre- 
vious fMRI data showing different responses in brain regions 
involved in the attachment and in the integration of auto- 
nomic states with social behavior to infant negative emotions 
between the mothers and non-mothers (e.g., hypothalamic- 
midbrain-limbic-paralimbic-cortical circuits) (Seifritz et al., 
2003; Feldman, 2007; Swain et al, 2008). 

The present study crucially extends the results obtained by 
our previous work (Ebisch et al, 2012) where, using the same 
paradigm, we provided consistent evidence for a synchrony in 
emotional response of mothers watching their own distressed 
child. Following that evidence, the present study tested in a larger 
sample the hypothesis that the above synchrony was modulated 
by the bond between the woman and the child. Therefore, a con- 
trol group of "other women" was recruited and compared with 
the "mother" group. A key improvement, from a methodologi- 
cal point of view, in the present work is represented by the use 
of a novel technique for extracting the thermal signal in pre- 
defined region of interest. The tracking algorithm allowed for 
an extraction of the signals that, despite the movements of the 
subjects, was objective over the entire experimental session, thus 
providing a stable and more accurate evaluation of tempera- 
ture variations, together with an increased number of processed 



frames per condition. The extraordinary potentialities of thermal 
infrared imaging in this field of research have been further con- 
firmed by this study, as we were able to record physiological 
measures of the vicarious reactions without interfering with the 
subjects' spontaneous behavior. 

Our study present some limitations that should be addressed. 
Firstly, the sample size is somewhat small. Therefore, our results 
should be considered preliminary, but suggestive of new insights 
in the study of the body communication between adults and 
children. Larger sample would be desirable, even though it has 
to be recognized that is it difficult, and somehow complex, to 
involve young children with their parents in a physiological study 
on children's distress. However, we plan to enlarge the sample 
size. The choice of preserving the ecological context determined 
the needs for the rejection of a not irrelevant number of dyads 
because of the excessive movements of the child that caused 
a artifacts that our tracking software could not fix. In addi- 
tion, it determined the impossibility for a obtaining a group 
average signal. 

To conclude, the findings show that the child's distress evo- 
cates in the observing women a spontaneous autonomic response, 
which could reflect an emotional sharing ability, stronger and 
faster in the mothers. Therefore, our study supports the hypoth- 
esis that the maternal bond can modulate the promptness of the 
adult's response to the child's needs. 
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